Fluctuation-dissipation theorem in an aging colloidal glass 
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We provide a direct experimental test of the fluctuation-dissipation theorem (FDT) in an aging 
colloidal glass. The use of combined active and passive microrheology allows us to independently 
measure both the correlation and response functions in this non-equilibrium situation. Contrary to 
previous reports, we find no deviations from the FDT over several decades in frequency (1 Hz-10 
kHz) and for all aging times. In addition, we find two distinct viscoelastic contributions in the aging 
glass, including a nearly elastic response at low frequencies that grows during aging. 

PACS numbers: 



Developing a statistical mechanical description of non- 
equilibrium systems such as glasses still remains an im- 
portant challenge in physics. One of the most interesting 
recent developments along these lines is the proposal to 
generalize the fluctuation dissipation theorem (FDT) to 
non-equilibrium situations [ij. The FDT relates the re- 
sponse of a system to a weak external perturbation to the 
relaxation of the spontaneous fluctuations in equilibrium 
0. The response function is proportional to the power 
spectral density of thermal fluctuations, with a prefactor 
given by the temperature. This suggests a generaliza- 
tion for systems out of equilibrium, in which the (non- 
equilibrium) fluctuations are related to the response via 
a time-scale-dependent effective temperature. While this 
has been studied extensively theoretically |3|], the exper- 
imental support for a meaningful effective temperature 
is unclear. There have been few experiments Q, 0, 
and the usefulness of the extension of the FDT to non- 
equilibrium situations is still a matter of controversy. 

Here, we introduce a combination of both active 
andpassive (fluctuation-based) microrheology techniques 
0, liot that provide a way to directly test the applica- 
bility of the FDT. We examine the validity of FDT in a 
colloidal glass, the synthetic clay of Laponite [ill H3. Il3| . 
For this system conflicting results have been reported 
previously 0, 0] , that may in part be due to the use of a 
limited experimental window in both frequency and ag- 
ing time. Here, we perform measurements over a wide 
range of frequencies and aging times. Contrary to previ- 
ous reports, we find no violation of the FDT and thus no 
support for an effective temperature different from the 
bath temperature. 

In addition, these measurements provide a new insight 
into the physics of the aging process. The microrheol- 
ogy done during the aging process allows us to explore 
the evolution of viscoelastic properties of the glass over 
a wider frequency range than hitherto explored, span- 



ning nearly 6 decades in frequency. The measurements 
reveal the existence of two distinct contributions to the 
viscoelasticity of the system: (i) a high-frequency vis- 
coelastic response in which the shear modulus increases 
rapidly with frequency; and (ii) a predominantly elastic 
(weakly frequency-dependent) response at lower frequen- 
cies, which becomes increasingly important as the system 
ages. 

The Einstein relation relates the diffusion of the parti- 
cle (i.e., position fluctuations) to its mobility. This is a 
special case of the FDT, and its generalization to a vis- 
coelastic system in Fourier space relates the power spec- 
tral density (PSD) of position fluctuations to the imagi- 
nary part of the complex response function a"(u>) 



|*H| 2 ) = ^«>) 



(i) 



Here, (|a;(u;)| 2 ) denotes the Fourier transform of the en- 
semble average (x(t)x(0)) . In a non-equilibrium system 
this suggests the introduction of an effective temperature 
in which T is replaced by T e g(u>) in Eq. lff|. 

To investigate the aging of colloidal glass, we study 
the motion of probe particles using optical tweezers. A 
custom-built inverted microscope |£| , equipped with two 
overlapping optical tweezers formed by two independent 
lasers (wavelengths A= 830 nm and 1064 nm ) focused 
to diffraction-limited spots. The latter drives the oscil- 
lations of the trapped particle with an Acousto-Optical 
Deflector, allowing us to measure the (active) response 
to a driving force. The (x, y) position of the particle is 
determined by a quadrant photo diode 0] with a spa- 
tial resolution of ~ 0.1 nm. The output signal from the 
photodiode is fed into a lock-in amplifier that measures 
the amplitude and phase of the particle displacement x(t) 
caused by an oscillatory motion of the drive laser focus. 
From the motion of the drive laser we determine the force 
F(t) acting on the particle. The response function is then 
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FIG. 1: The displacement power spectral densities (PSD) 
as a function for frequency of 1.16 fim silica probe particles 
with increasing age after preparing the sample. Fluctuations 
were recorded for 45 seconds with the 830 nm laser focus and 
results averaged in x-y directions for 1 bead several times. 
Aging times are given in the legend. The filled squares show 
the PSD of a bead in pure water for comparison. An acoustic 
noise signal around / ~ 200Hz is cut out from the curve at 
the latest stage of aging where the displacement signal was 
the lowest. All experiments were done at 21°C. 



given by a(u) = x u /F u , where and denote the 
fourier transforms of x(t) and F(t). By measuring the 
PSD of the same beads under the same conditions in wa- 
ter, we are able to calibrate both trap stiffness 0, 3| and 
particle displacement for active and passive data 8j . 

For the passive measurements the shutter in front of 
the driving laser is closed and the spontaneous fluctua- 
tions of the particle position are recorded for a minimum 
time of 45 s. From the displacement time series, we cal- 
culate the displacement power spectral density by Fast 
Fourier Transform ||. Comparing the response function 
from the active microrheology with the fluctuation spec- 
tra, we can directly check the validity of the FDT during 
the aging of the system, as well as resolve the frequency- 
dependent viscoelastic properties during the aging of the 
glass. 

The colloidal glass under study is a suspension of 
Laponite XLG in ultra pure water. After mixing the col- 
loidal particles with the water, the system spontaneously 
evolves from an initially liquid and ergodic state to a non- 
ergodic glassy state that exhibits elastic behavior [l2j |. 
For a particle concentration of 2.8 wt%, the rate of ag- 
ing is slow enough, on the one hand, that no significant 
structural and dynamic changes occur during each in- 
dividual active and passive microrheology measurement 
lasting at most 2 min. Nevertheless, the system evolves 
fast enough to allow us to follow the evolution from 'liq- 
uid' to 'solid' (the glass no longer flows when the sample 



cell is tilted) within about 8 hours. The dispersions are 
filtered (Millipore Millex AA 0.8 /an filter) to obtain a 
reproducible initial state This defines the zero of 

aging time t a = 0. Immediately after filtration, a small 
fraction (< 10~ 4 vol%) of silica probe beads (diameter 
1.16 /im ±5%) are mixed with the Laponite dispersion. 
The solution is then introduced into a sample chamber 
of about 50 (A volume, consisting of a coverslip and a mi- 
croscope slide separated by a spacer of thickness 70 //m. 
This is sealed with vacuum grease to avoid evaporation 
of sample. We then trap a single silica bead and perform 
the active and passive experiments on it. 

Since the system evolves toward a non-ergodic state, 
the time average may not be equal to ensemble average 
for the measured PSDs. To investigate this, we confirmed 
that reproducible PSDs were obtained for the same aging 
time, independent of bead position, during all stages of 
aging. We also confirmed that our results do not depend 
on the time interval used to compute the time average. 
Thus, we can use the time-averaged PSD without aver- 
aging over several beads in our study. Figure ^ shows 
the (passive) displacement PSD for different aging times. 
It is evident that the particle motion progressively slows 
down with increasing aging time, reflecting the increase 
in viscosity of the system. Qualitatively two regimes of 
aging are seen: for t a < 200 min the PSD can be de- 
scribed by a single power law. At longer aging times two 
distinct slopes appear in the log- log plot (Fig. [TJ. 

We also measure the (active) response of the same 
bead used in passive measurements, as a function of 
aging time and for oscillation frequencies of / =1.2, 
10.8, 116, 1035 and 12000 Hz. To directly compare the 
(passive) fluctuations with the (active) response, we ex- 
press our fluctuation PSDs normalized in such a way 
as to permit a direct comparison with the measured a 
in the form of Eq. (JTJ- Thus, we plot the measured 
PSD multiplied by ujj (2&bT). We obtain the real part 
using a Kramers-Kroning (principal-value) integral 



d£. The cutoff error due to a fi- 



nite range of frequencies sets an upper limit to a'(tu) 
about a decade lower than that of a"(uj). Fig. |21 depicts 
the real and imaginary parts obtained from the active 
and passive methods at an early (t a = 100 min) and a 
late stage of aging (t a = 300 min). We see that the re- 
sults are identical to within the experimental accuracy, 
showing that there are no deviations from the FDT in 
this system over the range of frequencies and aging times 
probed in our experiments. Note that the small devia- 
tions between the respective a'(u>) values at high frequen- 
cies are likely due to cut-off errors in the Kramers-Kronig 
integrals because the drawn lines, obtained by extrapo- 
lating the measured power law of a" to infinity, show 
significantly better agreement. Since we have directly 
compared both the real and imaginary parts of the re- 
sponse functions, this represents a stronger test of the 



3 




k> (rad/s) co (rad/s) 

FIG. 2: Comparison of active and passive results: Real a!(u)) 
and imaginary a"(ui) at t a = 100 and 300 min obtained from 
active (filled symbols) and passive (open symbols) microrhe- 
ology performed on the same 1.16 /im diameter silica bead 
in the same sample. For the passive experiments, the imagi- 
nary parts of the response functions are obtained directly and 
real parts are calculated with a Kramers-Kronig integral. The 
lines show the fits to Eq. @ . At early stages of aging the data 
can be described with one power law, while at later stages, 
a superposition of two power laws is needed to describe the 
whole frequency range. The amplitude of oscillation for the 
active experiments was 77 nm. 



Teff/Tbath 


ta 


7.5 rad/s 


68 rad/s 


728 rad/s 


6.5 krad/s 


75 krad/t 


0-2 h 
2-4 h 
4-6 h 
6-8 h 


0.75 ± 0.3 
1.2 ± 0.3 
1.4 ± 0.3 


1± 0.1 
1 ±0.1 
1 ±0.1 
0.85 ±0.1 


0.95 ± 0.1 
1 ± 0.1 
1.1 ± 0.1 
1.0 ± 0.1 


0.85± 0.1 
0.9 ± 0.1 
1.1± 0.1 
1.1 ± 0.1 


1.0 ± 0.1 

1.0 ± 0.1 

1.1 ± 0.1 
1.0 ± 0.1 



TABLE I: The effective temperature obtained for different 
frequencies averaged over 2 h time intervals. Within the un- 
certainty in the experiments, T^/T^^ = 1 



FDT than previous measurements |(J and demonstrates 
that the FDT holds in this non-equilibrium system. 

In Fig. 01 we plot the extracted a" as a function of ag- 
ing time for several different frequencies. As can be seen 
the active and passive data agree very well. This figure 
confirms again that to within the experimental uncer- 
tainty the FDT holds: the measured effective tempera- 
ture does not differ from the bath temperature. The re- 
sulting effective temperature T cff /T heith = ap' assivc /a^, tive 
is shown in the Table I. We conclude that the measure- 
ments show the FDT in the form of the Einstein relation 
is valid for all frequencies probed here, and can be used 
for all the stages of aging in this system. The method 
also allows us to obtain the viscoelastic properties over a 
very wide frequency range; classical (macroscopic) rheol- 




FIG. 3: The comparison of a" extracted from passive (open 
symbols) and active (filled symbols) measurements as a func- 
tion of aging time. For the lowest frequency, we could not 
measure reliable data longer than 300 min, since the signal to 
noise to ratio became of the order of 1 with increasing aging 
time, as the material became stiff. 

ogy is limited to frequencies up to about 100 rad/s [l2| . 
In Fig. we observe a gradual decrease of the PSD for 
higher frequencies and a more rapid change at lower fre- 
quencies. The response function is directly proportional 
to the PSD which in turn should be inversely related to 
the complex shear modulus G* = G'+iG". With increas- 
ing aging time, as the sample becomes more viscous and 
solid-like, one would expect that both the elastic modulus 
G' and the viscous modulus G" increase. This is consis- 
tent with the decrease of the PSD and consequently of 
the response function. As was mentioned above, at late 
stages of aging two distinct slopes appear in the PSD 
(Fig. 2]) . This suggests the existence of two distinct con- 
tributions to the viscoelasticity during aging. Assum- 
ing the generalized Stokes formula for the viscoelastic 
response function 0, @ > we obtain excellent fits to the 
data assuming a simple addition of two power law con- 
tributions to the complex shear modulus (only a single 
power-law contribution at the early stages of aging) : 



v ' QttRG*{uj) Ci(-iuj) a + C 2 {-iLu) b y ' 

The fit of the imaginary part of the response function 
from the passive measurements with the imaginary part 
of the above functional form is shown by the drawn lines 
in Fig. [21 To demonstrate the quality of this model for 
describing the data, we also plot the real part with the fit- 
ting parameters obtained from the imaginary part. The 
agreement is remarkable, especially with the active data 
at high frequencies. Fig. 0] depicts the evolution of the 
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FIG. 4: For t a < 150 min, the complex response function 
can be described by single power law of the form l/(— iuj) a . 
After about 150 min, a slowly relaxing contribution emerges 
leading to a response function of the from in Eq. (2). 



fitting parameters, i.e., the exponents of power laws and 
the weight factors for the contribution of the two vis- 
coelastic contributions as a function of aging time. The 
exponent and amplitude of one of the components do not 
change with aging time while the amplitude of the other 
one grows appreciably for aging times longer than about 
250 min. 

These results demonstrate the existence of two distinct 
contributions in the viscoelasticity of the system. In ad- 
dition to a strongly frequency-dependent viscoelastic re- 
sponse at high frequencies, we also observe the slow devel- 
opment of a more elastic (weakly frequency-dependent) 
response during the aging. In fact, this appears to be the 
main characteristic of the aging in this system. A simi- 
lar description in terms of a network in a more fluid-like 
background has been suggested before for polymeric gels 
our case, this can be attributed to the growth 
of a tenuous network-like structure, in addition to the 
more viscous response that is always present. We note, 
however, that although our system also becomes non- 
ergodic there is no evidence for spatial inhomogeneity, 
neither in our system, nor in glassy systems in general. 
Recently, it has been suggested that glassy systems may 
be dynamically heterogeneous an d that one should 
look into the correlation of spatial and temporal dynam- 
ics in order to detect the heterogeneity. Therefore the 
homogeneity of the combined spatial and temporal dy- 
namics merit further detailed study. 

In summary, we see a good quantitative agreement be- 
tween the response function and the spontaneous thermal 
fluctuations, implying that we observe no violation of the 
FDT in this non-equilibrium system. Equivalently, we 
find an effective temperature that does not differ from 
the system temperature. It is important to note that 



these measurements provide a much more direct test of 
the FDT than prior experiments, since we directly mea- 
sure the response and the corresponding fluctuations over 
the same and very large range of frequencies. 
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